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We investigate quantum charge and spin pumping in armchair graphene nanoribbons under a sin- 
gle ac gate voltage connected with nonmagnetic/ferromagnetic leads via the nonequilibrium Green's 
function method. In the case of nonmagnetic leads, where only part of the nanoribbon is subject 
to an ac gate voltage to break the left-right spatial symmetry, we discover that peaks of the charge 
pumping current appear at the Fermi energies around the subband edges in the ac-field-free region 
of the nanoribbon. In the case of ferromagnetic leads with the lead magnetizations being antiparallel 
to break the left-right symmetry, similar peaks appear in the spin pumping current when the Fermi 
energies are around the edges of the the majority-spin subbands in the ferromagnetic leads. All 
these peaks originate from the pronounced symmetry breaking in the transmissions with energies 
being around the corresponding subband edges. Moreover, we predict a pure spin current in the 
case of ferromagnetic leads with the whole graphene nanoribbon under an ac gate voltage. The 
ac-field-strength and -frequency dependences of the pumping current are also investigated with the 
underlying physics revealed. 
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I. INTRODUCTION 

Graphene and its lower dimensional cousins, graphene 
nanoribbons (GNRs) and carbon nanotubes, exhibit 
abundant new physics and potential applications and 
hence have attracted much interest in recent years ji"— 
Among different works in this field, the effect of an ac 
field on the electrical and optical properties in these ma- 
terials is one of the main focuses of attention. Many 
interesting phenomena have been reported, such as 
the photon- assisted transport,^' the photovoltaic Hall 
effect,— the dynamical Franz-Keldysh effect^^ and the 
quantum pumpingj^^"— 

The quantum pumping, which is highly related to the 
ratchet effeclj^"— and the photogalvanic effect,—"— de- 
scribes the generation of a direct current at zero bias 
in a spatially asymmetric system under ac fieldsi^>2^ 
The quantum pumping of charge and spin currents has 
been investigated theoretically^"— and observed experi- 



mentally in semiconductor quantum dots 
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quantum 



wirea^"— and also in carbon nanotubesJ^ Recently, the 
quantum charge and spin pumping in GNRs has also 
aroused growing attention. ^^"^^ Most of these studies 
focus on the pumping involving more than one time- 
dependent parameters,-W:"-2 partially because the single- 
parameter pumping is only possible beyond the adiabatic 
approximation and thus needs more complex theoreti- 
cal toolj^ Nevertheless, the single-parameter pumping 
is more favorable to the application, since the reduction 
in number of necessary contacts makes the scalable and 
low-dissipative device more promisingjS 

So far, there are few investigations on the single- 
parameter pumping in the GNR. Torres et al.^^ studied 
the charge pumping in the ribbon of small size, where 



the behaviour is dominated by the resonant tunneling. 
However, the single-parameter pumping in a typical one- 
dimensional GNR, i.e., with large length and small width, 
has not yet been investigated. Moreover, to the best of 
our knowledge, there is no work on single-parameter spin 
pumping in this system. The aim of this study is to fill 
these spaces. 

In this paper, we present a detailed study of the single- 
parameter quantum charge and spin pumping in arm- 
chair GNRs contacted with nonmagnetic/ferromagnetic 
leads via the nonequilibrium Green's function method^ 
We first address the case of nonmagnetic leads, in which 
only part of the GNR is subject to an ac gate voltage 
to break the left-right spatial symmetry. It is discovered 
that peaks of the negative (positive) charge pumping cur- 
rent appear at the Fermi energies around the energy max- 
imums (minimums) of the subbands in the ac-field-free 
region of the GNR. Then we turn to the case of ferromag- 
netic leads with the lead magnetizations being antiparal- 
lel in order to break the left-right symmetry. It is shown 
that a pure spin current can be achieved when the whole 
GNR is under an ac gate voltage. We also predict peaks 
in the negative spin pumping current at the Fermi ener- 
gies around the energy maximums of the majority-spin 
subbands in the ferromagnetic leads. In the appendix, we 
discuss the time-dependent ballistic magnetotransport in 
armchair GNRs contacted with ferromagnetic electrodes 
and show that the results under the cutoff-energy ap- 
proximation, artificially introduced by Ding et al.r^ are 
qualitatively different from those obtained from the exact 
calculations. 

This paper is organized as follows. In Sec. II, we set 
up the tight-binding Hamiltonian and the formula of the 
pumping current. The numerical results are presented in 
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Sec. III. Finally, we summarize in Sec. IV. 



II. MODEL AND FORMALISM 



We consider an armchair GNR with a single gate 
voltage applied between two nonmagnetic/ferromagnetic 
leads as shown in Fig.[T] In this system, the tight-binding 
Hamiltonian can be written as 



H — Hl + Hr + Hg + Ht , 



(1) 



in which Hl {Hb) represents the Hamiltonian of the left 
(right) lead, Hg stands for the Hamiltonian of the GNR 
and Ht describes the hoping between the GNR and the 
leads. These terms can be written as 



Hl = 
Hb = 
+ 

= 



c], 



E ^Lac\^^Cj^„, (2) 



iR.cr 

^M_Rsin6'c|^^Q^_^- E ^«'^4r<7Cjh<t,(3) 
E^'«<-^'«-- E hdi.d.^o, (4) 

-EE ^Ta^aCjc.+H.C. (5) 

a—L,R {ig ) ,(7 



Here Ci^^ (digcr) and cj^^. (c?J g.) are the annihilation and 
creation operators of the electron with spin a on lattice 
site ia (ig) in the leads (GNR); {i,j) denotes pair of near- 
est neighbors; Eaa = Eao + crMa are the on-site energy 
of the spin-up (cr = -I-) or spin-down (cr = — ) band in the 
ferromagnetic leads with Ea+ = E^- = Eao for the non- 
magnetic leads; 9 is the angle between the magnetization 
directions of right and left leads; taa and txa- represent 
the hoping parameters in the leads and between the leads 
and the GNR, respectively, which are set to be equal to 
the hoping parameter tg = 2.7 eV in the GNR, unless 
otherwise specified. The on-site energy of the GNR Vi^ 
takes the values of -I- VacCOs{Qt) and Vg in the ac- 
field-applied and -free regions of the GNR, respectively, 
with Vac and fl being the magnitude and frequency of the 
applied ac gate voltage. Here we apply two static gate 
voltages Vg and Vg to modulate the on-site energies in 
these two regions of the GNR independently to facilitate 
the identification of the influence of the pumping from 
different regions. 

Exploiting the nonequilibrium Green's function 
method^ the time-averaged current can be written as 
(the current flowing from left to right is deflned to be 
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FIG. 1: Schematic view of the armchair GNR connected with 
two nonmagnetic/ferromagnetic leads. 



positive) 

I = 



,i{s' —s)t 



e r" , f°° dede'dei 
/ at L 

X Tr{G''(£,£i)ffl(ei)G'-(£i,e')rL(e')/L(e') 
- G-{e,e^)tL{ei)G^ei,e')rB.{e')fB{e')} . (6) 

Here Tq = 27r/f}; ^^("^^(e, e') is the retarded (advanced, 
lesser) Green's function in the ac-fleld-applied region 
of the GNR; /a(e) is the Fermi distribution; Tai^) = 
i[S;(£) - SJi(e)], with I]2^°^(£) and T.''^''\e) representing 
the retarded (advanced) self-energies from the left lead 
and both the right lead and the ac-field-free region of the 
GNR, respectively. The self-energy has the form 



s°+ 
s5_ ' 



(7) 



where □ = r, a, the rotational matrix Rq, is defined as 
f cos % 



Ra — 



sm ■ 



V sin ■ 



(8) 



with eL=0 and 0^ = 61. 

Then we need to calculate G"'(£,£'). It is known that 
the eigenstates of the isolated GNR with and without an 
ac gate voltage, labelled as $(i) and $o(^) respectively, 
satisfy the relatio n^^i^" 



$(i) = $o(Oexp|--y KcCos(m)dT|. (9) 

Thus the Green's Function of the isolated ac-field-applied 
GNR (£,£') has the form 



g'^{sr + nfl, e[. + m£L) = 2nS{er — £j,).g'"(£r, n, m), (10) 

Vac \ J j Vac 



g^'iSr, n, m) = ^ Jn-N{^)Jm-N{^)go{£r + Nil). 



N 



(11) 



In these equations Er € [—§,§), 5o(^) is the corre- 
sponding Green's function of this GNR without the ac 
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gate voltage, which can be obtained by the recursive 
method;^ From the Dyson equation 

G'-(e,£')= 5'' (£,£') + / EOS'- (£i)G'- (£!,£'), 

(12) 

one obtains 

G"'(£r + nn, e; + mQ) = 27r5(e^ - e'^)G''{er, n, m), (13) 
where G {£r,n,m) is determined by 

CT'{er,n,m) = g^jsr, n, m) + g''(£r, n, ni) 

ni 

y.t''{er + nin)(f{er,ni,m). (14) 
Substituting Eq. into Eq. ([5]), / can be written as 

/oo 
rf£[rrfl.(£)/L(£) - n^MfR{e)l (15) 

in which 

T2rAs) = Tr{fL,(£)G°(e,e + nf])ffl,(e + nr!) 

X G^(e + nf7,e)} (16) 

is the transmission probability from left to right leads 
involving the absorption {n > 0) or emission {n < 0) of 
\n\ photons of electrons with initial energy e and spin a. 
Here we have limited ourselves in the parallel and an- 
tiparallel configurations of the electrode magnetizations, 
and thus the contribution from different spin bands can 
be calculated separately. 

In the spatially asymmetric system, generally speaking 
^Lfia-(^) 7^ ^^RLai^)- Thus cvcn whcn the external bias 
is absent, i.e., /l(£) = fni^), the time-average current 
can be nonzero and a pump current emerges. In the zero 
temperature limit, the pumping current with spin a is 
given by 

^pump -tT. r '^^^'^LRAe) - Ti^Laie)], (17) 
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FIG. 2: (Color online) Armchair GNR connected with non- 
magnetic leads. Nf, — 100. (a) Charge pumping current 
as function of the Fermi energy Ep for Vg = Vg — (red 
solid curve), Vg = 0.05tg, Vg = (blue dashed curve) and 
Vg — 0, Vg — 0.05tg (green dotted curve). Note that the 
scale of Ef for Vg = and Vg — 0.05tg is on the top of the 
frame. The vertical black dashed lines indicate the energy 
maximums and minimums of the subbands in the ac-field- 
free region of the GNR. (b) The schematic illustration of one 
peak in the negative pumping current for Vg = Vg = 0. The 
black dashed line indicates the Fermi energy. The red and 
blue solid arrows represent the transmissions T}^b.ct{^) s-nd 
TftLai^) in the energy regime satisfying both Ep — fl < e < Ep 
and iJniax ~ ^ < e < £max- The red and blue dashed ar- 
rows represent the transmissions Tlj^^le) and T^l^{e + 
for e < Ep — Q,, which cancel each other. 



with E-p denotes the Fermi energy. The charge and spin 
pumping currents are defined as /pump = -^ilump + -^pump 
and /^ump = -^i^ump - -^iTump. respectively. Further con- 
sidering the time- reversal symmetry, one has T£j^^(e) — 
T^^^{e + nil). This means that the photon-assisted 
transmission with the initial and final energies both be- 
low the Fermi energy is cancelled by the corresponding 
one in the opposite direction and hence cannot contribute 
to the pump current. Consequently, Ip^^p can be ex- 
pressed as 

Cmp = r E / de[T£j,^{e) - T^^Ae)]- (18) 



III. NUMERICAL RESULTS 



In this section, we present the numerical results of the 
quantum charge and spin pumping currents in the arm- 
chair GNR connected with nonmagnetic/ferromagnetic 
leads. In our computation, the parameters are chosen 
to be Nw = 41, Nac = 400, V^c = ^ ^ O.Olig and 
Vg = Vg = unless otherwise specified. Since the width 
satisfies Nw — 3M-I-2 with M being an integer number, 
thus this armchair GNR is metallici^ Also note that the 
length of the GNR is large enough to totally suppress the 
contribution from the evanescent modesi^ 
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A. Quantum pumping in armchair GNR connected 
with nonmagnetic leads 

We first study the quantum charge pumping in the 
armchair GNR connected with nonmagnetic leads, where 
the on-site energies are chosen to be = E^- = 0. 
We set a finite length of the ac-field-free region of the 
GNR, i.e., Nb = 100, in order to break the symmetry be- 
tween the transmissions T£j^^{e) and ^^'^ hence 
induce the pumping current.--'- - In Fig.[5fa), we plot the 
charge pumping current /pump against the Fermi energy 
Ep for Vg =Vg = (red solid curve). Here and hereafter, 
the energy zero point is chosen to be the Dirac point 
of the pristine GNR. It is shown that peaks of negative 
(positive) current appear around the energy maximums 
(minimums) of the subbands in the GNR indicated by 
the vertical black dashed lines. 

These peaks are understood to originate from the pro- 
nounced symmetry breaking in the transmissions with 
energies around the subband edges in the ac-field-free 
region. We take one peak in the negative current as 
an example to illustrate this physics in Fig. [2][b) . Here 
we present the relevant subbands in the ac-field-applied 
and -free regions of the GNR, whose energy maximums 
are labelled by i?max and E'^^^, respectively. The rele- 
vant subbands in the left and right leads are also shown, 
in which propagating modes exist in the whole energy 
range investigated here. The transmissions T^j^^{e) and 
Tj^Raie) with £^max — < £ < -^max are plotted in this 
figure as blue and red solid arrows, respectively. It is seen 
that the initial (final) energies of these transmissions are 
lower (higher) than E[-^^^. Although both transmissions 
are allowed in the ac-field-applied region due to the side- 
band effect, Tlj^^{e) is forbidden in the subband in the 
ac-field-free region owing to the lack of the propagat- 
ing modes above E^^^. This makes the asymmetry in 
these transmissions become notable. It is also noted that 
only for Ep — fi < e < Ep, the single-photon assisted 
transmissions can contribute to the pumping current as 
indicated by Eq. ((T5)) . As a result, a large negative pump- 
ing current appears when iS^iax — < Ep < E'max + ^ 
and reaches its peak value at the Fermi energy around 
^max- Similar symmetry breaking can also be found in 
the multi-photon assisted transmissions. After taking 
these transmissions into account, the peak of the pump- 
ing current appears at almost the same Fermi energy but 
with wider width. The presence of the peaks in the pos- 
itive pumping current is based on the similar physics. 
When the Fermi energy is around the minimum of one 
subband in the ac-field-free region, this subband only 
contributes to the transmissions from left to right, but 
not to the opposite ones. This symmetry breaking leads 
to the peak of positive current around this minimum. 

In order to further verify the above physics, we also 
plot /pump with the modulated on-site energy Vg = 
0.05tg, v[ = (blue dashed curve) and Vg = 0, V^ ^ 
0.05tg (green dotted curve) in Fig. [5Ka). It is seen that 
in the former case, the Fermi energies corresponding to 
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FIG. 3: (Color online) Armchair GNR connected with non- 
magnetic leads. Nb = 100. (a) Charge current as function of 
the ac-field strength Vac with the frequency Q — O.Oltg for 
Ef — 0.065tg (red curves with •) and 0.132tg (blue curves 
with □). The solid and dashed curves represent the cur- 
rents with no bias (i.e., the pumping current) and negative 
bias = Ef - and E^ = Ef + fl/l6. The black 

dashed line indicates the position of zero current, (b) Charge 
pumping current as function of the ac-field frequency Q with 
Vac ~ O.Oltg for Ef — 0.065tg (red solid curve with •) and 
0.132tg (blue dashed curve with □). 



peaks of negative current are the same as those in the 
unmodulated case; whereas in the latter one, the Fermi 
energies corresponding to peaks of positive current are 
higher than the unmodulated ones by 0.05ig. These be- 
haviours agree with the physics discussed above: it is the 
subband edges in the ac-field-free region that determine 
the energy range where the pronounced symmetry break- 
ing appears, and hence the Fermi energies corresponding 
to the peaks in the pumping current. Moreover, the mod- 
ulation of the on-site energy also brings about more in- 
teresting phenomena. It is shown that in the case with 
Vg = Omtg and Vg' = (Vg = and V^ = 0.05ig), the 
peaks of negative (positive) current appear with larger 
magnitude than those in the unmodulated case, while all 
peaks of positive (negative) current disappear. These be- 
haviours for Vg — O.Q5tg and Vg = can be understood 
as follows. When the Fermi energy is around i?max shown 
in Fig.[2I^b), the relevant transmissions from right to left 
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in this case are stronger than those without modulation, 
since the increase of i?max makes more states in the sub- 
band in the ac-field-apphed region contribute to these 
transmissions through the sideband effect. Consequently, 
the asymmetry in the transmissions becomes more pro- 
nounced and thus the magnitude of the peak in the neg- 
ative current becomes larger. On the other hand, when 
the Fermi energy is around the minimum of one subband 
in the ac-field-free region, the relevant transmissions in 
both directions are forbidden in the corresponding sub- 
band in the ac-field-applied region, because the energy 
minimum of this subband is much higher than the related 
energies of transmissions. This cancels the asymmetry in 
the transmissions and hence suppresses the peak in the 
positive current. The behaviours in the case with Vg — 
and Vg = 0.05ig can be understood in the similar way. 

We stress that the above picture is only valid in the 
long ribbons, where the energy spectrum can be consid- 
ered quasi-continuous, but invalid in the short ribbons. 
In the latter case, the field frequency is smaller than the 
energy level spacing, thus the behaviour of the pumping 
current is dominated by the resonant tunnelingfS£ instead 
of the physics presented above. 

Then we turn to the ac-field-strength and frequency 
dependences of the pumping current. In Fig. ^a), the 
pumping current is plotted against the ac-field strength 
Vac with frequency ft — 0.0 It^ for two Fermi energies 
Ep = O.OGbtg and 0.132tg, which correspond to the two 
peaks of positive current in Fig. [Ifa). It is seen that 
-^pump increases monotonically with Vac for various Fermi 
energies. This is because the photon-assisted transmis- 
sions are enhanced when the ac field becomes stronger. 
Besides the pumping currents, i.e., the current with no 
bias, the currents with negative bias = Ep — U/16 and 
Ep = Ep + n/16 are also plotted in this figure by dashed 
curves. It is seen that the current is negative in the ab- 
sence of an ac field. When the ac field is strong enough, 
the pumping current can surpass the current driven by 
the negative bias and hence make the net current become 
positive, i.e., along the direction opposite to the external 
bias. 

In Fig. ^h), we plot the pumping currents against 
the field frequency with strength Vac = O.Oltg for 
Fermi energies Ep — 0.065tg and 0.132tg. It is seen 
that /pump presents a nonmonotonic frequency depen- 
dence when n < Vac- The scenario is as follows. On 
one hand, with an increase of the frequency, Vac/^ de- 
creases and thus the photon-assisted transmissions be- 
come weaker. On the other hand, the energy range of 
the transmissions contributing to the pumping current 
increases with as shown in Eq. (|18p . The competition 
of these two factors leads to the complex frequency de- 
pendence. For r2 > Vac, the pumping current is shown 
to decrease with increasing frequency. This is because 
the ac field becomes ineffective in exciting photons when 
frequency is high enough. 
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FIG. 4: (Color online) Armchair GNR connected with ferro- 
magnetic leads with the lead magnetizations being antiparal- 
lel. (a) Spin pumping current /pump as function of the Fermi 
energy Ef for El+ = 0, Er+ — —3tg (red solid curve) and 
—S.ltg (blue dashed curve). A'^i, = 0. Note that the scale 
of Ef for Er+ = —3.1tg is on the top of the frame. In the 
inset, the spin-up subbands in the GNR and the ferromag- 
netic leads relevant to the peaks of the pumping current at 
Ef = (red curves) and 0.126tg (blue curves) are schemat- 
ically plotted. The solid curves represent the subbands in 
the case with El+ = and Er+ = —3tg, while the dashed 
ones in the right panel represent the subbands in the right 
lead for /5h+ = —S.ltg. (b) Charge and spin pumping cur- 
rents together with the spin-up and -down pumping currents 
against the Fermi energy Ef for El+ ~ and Er+ = —3tg. 
Nb = 100. In these figures, the vertical black dashed (pink 
dotted) lines indicate the energy edges of the subbands in the 
ac-field-free region of the GNR (the majority-spin subbands 
in the ferromagnetic leads). 



B. Quantum pumping in armchair GNR connected 
with ferromagnetic leads 

Now we investigate the quantum charge and spin 
pumping in the armchair GNR connected with ferromag- 
netic leads. In order to break the left-right symmetry, the 
magnetizations of the two leads are set to be antiparallel, 
i.e., ELa = Eji-rr, and hence the ac-field-free region of 
the GNR becomes unessential. We first focus on the case 
without this region, i.e., Nh = 0. In this case, the system 
satisfies the symmetry Rt^TZlr, where /?^ is the spin ro- 
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tation operator defined by Eq. ([8]) with rotation angle tt 
and TZlr is the spatial reflection operator in the direction 
from left to right. Thus, the spin-up and -down pump- 
ing currents always satisfy the relation /pJu,np — —^pump- 
Consequently, the charge pumping current vanishes and 
a pure spin current can be achieved. 

We plot the spin pumping current /pu^p fo^' -E'L+ = 
and = — 3tg as red solid curve in Fig. IDJa). It 

is shown that peaks of negative current appear at the 
Fermi energies around the maximums of the majority- 
spin (L— and R+) subbands in the ferromagnetic leads, 
as indicated by the vertical pink dotted lines. In order to 
reveal the underlying physics, we schematically plot the 
spin-up subbands in the GNR and the ferromagnetic in 
the inset of Fig. lUJa). Here only the subbands relevant 
to the peaks at iSp = (red sohd curve) and 0.126tg 
(blue solid curve) are shown. Comparing this inset with 
Fig. m^b) , one reaches the conclusion that the physics of 
these peaks is similar to that in the case of nonmagnetic 
leads: they come from the pronounced asymmetry in the 
transmissions with energies around the maximums of the 
majority-spin subbands in the ferromagnetic leads. In 
order to further verify this physics, we low the on-site 
energy Eji+ by O.lig and plot the corresponding /pump 
as blue dashed curve in Fig. Hl^a). One observes that 
the Fermi energies corresponding to peaks also decrease 
by O.ltg as expected. The only exception is that the 
peak previously at Ep = 0.126tg is absent in this case. 
This feature is based on the similar physics leading to 
the absence of the corresponding peaks in the case of 
nonmagnetic leads. As shown in the inset of Fig.|4l[a), the 
energy gap between the subbands relevant to the absent 
peak in the GNR (blue solid curve in the middle panel) 
and right lead (blue dashed curve in the right panel) is 
much larger than the photon energy. Consequently, these 
subbands cannot contribute to the relevant transmissions 
in either direction. This leads to the absence of the peak. 

Now wc turn to the case with finite N),. Due to the 
breaking of the symmetry R.^TZlr, both charge and spin 
pumping currents can be finite. In Fig. IHb), we plot 
these two pumping currents against the Fermi energy for 
Elj^ = and Erj^ — —itg as red solid and blue dashed 
curves, respectively. By comparing these results with 
those in Figs.dja) andlD^a), similar features can be seen 
between the spin or charge pumping currents here and in 
the case of ferromagnetic leads for A^{, = or the case of 
nonmagnetic leads, i.e., peaks in the spin (charge) pump- 
ing current appear at the Fermi energies around the edges 
of the majority-spin subbands in the leads (the subbands 
in the ac-field-free region of the GNR) . The scenario is as 
follows. Here the pumping current comes from two kinds 
of symmetry breaking. The first is due to the antiparallel 
magnetizations of the ferromagnetic leads, which induces 
peaks where the Fermi energies are around the edges 
of the majority-spin subbands in the leads, as discussed 
above. From Fig. IH^b), one further observes that these 
peaks have similar magnitude but the opposite sign in 
the spin-up (yellow chain curve) and -down (green dot- 



ted curve) pumping currents. This results in peaks in 
the spin pumping current but provides negligible contri- 
butions in the charge pumping current. The second kind 
of symmetry breaking is from the presence of the ac-field- 
free region of the GNR, which induces peaks where the 
Fermi energies are around the subband edges in this re- 
gion, as demonstrated in Sec. IIIA. These peaks tend to 
cancel each other in the spin pumping current, but man- 
ifest themselves in the charge pumping current. Under 
the joint effect of these two kinds of symmetry breaking, 
the behaviour of the spin or charge pumping current re- 
sembles that in the case of ferromagnetic leads for A^^ = 
or the case of nonmagnetic leads. Moreover, some new 
features are induced due to the interplay of these two 
kinds of symmetry breaking, such as the peaks in the 
negative spin pumping current around E-p = 0.075ig and 
the peak in the positive charge pumping current around 
E-p — Q.2A7tg. In addition, we verify the ac-field-strength 
and frequency dependences of the spin and charge pump- 
ing currents with ferromagnetic leads resemble those of 
the charge pumping current with nonmagnetic leads and 
hence do not repeat here. 



IV. SUMMARY AND DISCUSSION 

In summary, we have performed a detailed study of 
the single-parameter quantum charge and spin pump- 
ing in the armchair GNR connected with nonmag- 
netic/ferromagnetic leads via the nonequilibrium Green's 
function method. We first study the charge pumping 
in the case of nonmagnetic leads. Here only part of 
the GNR is subject to an ac gate voltage in order to 
break the left-right spatial symmetry. We discover that 
peaks of the negative (positive) charge pumping current 
appear when the Fermi energies are around the energy 
maximums (minimums) of the subbands in the ac-field- 
free region of the GNR. This phenomenon comes from 
the pronounced symmetry breaking in the transmissions 
with energies being around these subband edges. We 
also discuss the pumping current with the modulated 
on-site energies in ac-field-free and -applied regions of 
the GNR. Our results show that the Fermi energies cor- 
responding to peaks are influenced by the energy-level 
shifting in the ac-field-free region, but independent of 
that in the ac-field-applied region. This is in agreement 
with the physics presented above. Furthermore, we find 
that the modulation of the on-site energies also induces 
the absence of some peaks. This is because the rele- 
vant transmissions are forbidden due to the lack of the 
propagating modes in the corresponding subbands in the 
ac-field-applied region. 

The ac-field-strength and -frequency dependences of 
the pumping current are also investigated. We show that 
the pumping current increases monotonically with the ac- 
field strength due to the enhanced photon-assisted trans- 
missions. Moreover, we discover that at low ac-field fre- 
quency, the pumping current exhibits a nonmonotonic 
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frequency dependence. This is due to the competition 
of the weakened photon-assisted transmissions and the 
increasing energy space of the states contributing to the 
pumping current. At high frequency, the pumping cur- 
rent decreases monotonicahy with increasing frequency, 
since the ac field becomes ineffective in exciting photons. 

More interesting features are seen in the charge and 
spin pumping in the case of ferromagnetic leads. Here we 
set the magnetizations of the two leads to be antiparallel 
to break the left-right symmetry. In the case with the 
whole GNR under an ac field, we show that the charge 
pumping current vanishes and a pure spin current can 
be achieved. We also find that peaks of negative spin 
pumping current appear at the Fermi energies around 
the maximums of the majority-spin subbands in the fer- 
romagnetic leads. The physics is similar to the pumping 
peaks with nonmagnetic leads. In the case with only part 
of the GNR subject to an ac gate voltage, both the charge 
and spin pumping currents can be finite. Our calculations 
show that the behaviour of the spin pumping current is 
similar to that in the case of ferromagnetic leads with the 
whole GNR under an ac field, while the behaviour of the 
charge pumping resemble that with nonmagnetic leads. 

Finally, we stress that the main features in the charge 
and spin pumping currents predicted here do not depend 
on the specific properties of the leads. Thus they ap- 
pear not only in the case of the simple-squared-lattice 
leads used in our model, but also in the case of the other 
types of leads, e.g., the two-dimensional superlattic e^ -i ' — ' 
and graphene leads. This provides more choices for the 
experimental investigations. 
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Appendix: Magnetotransport in graphene under ac 
gate voltage 

Recently, Ding et al^ studied the magnetotransport 
of a similar structure, i.e., graphene connected with fer- 
romagnetic leads in the presence of an ac gate voltage. 
Similar to our model, they treated the transport inside 
the graphene and leads as the ballistic transport, i.e., 
without considering any scattering. However, their treat- 
ment of the interface between graphene and leads is very 
different. They neglected the momentum dependence of 
the coupling matrix and replaced all coupling matrix el- 
ements by only one phenomenological parameter. This 
leads to a divergence in the transmission. To remove the 
divergence, they artificially introduced a cutoff energy. It 



is exactly due to this cutoff energy, many pronounced fea- 
tures in the time-dependent magnetotransport are pre- 
dicted in that paper. 

In the following, we will demonstrate that these pro- 
nounced features are just the artificial results induced 
by the inappropriate treatment of the interface and the 
introduction of the cutoff energy. Since there is no cut- 
off energy in our model, we can calculate the transport 
properties without any approximation used in Ref. Issl . 
and compare the obtained results with corresponding one 
with the cutoff energy. In the main text of this investiga- 
tion, the parameters of the ferromagnetic leads are set as 
El^ El- and tra = tg. In this appendix, in order to 
reveal the problems in that paper more clearly, we change 
the parameters to El+ = El- = and txa = tg^/T+lrP 
with P = 0.8 to make Vl+ITl- = (1 + P)/{1 - P) as 
done in Ref. [s^. In addition, we only discuss the case 
with the whole GNR under the ac gate voltage, i.e., 
Nb — 0. Thus for the total transmission Tlr/rl{£) = 
T^anT^R^/RLai^)' ^LRis) = TflL(e) duc to the Symme- 
try of this system. The other parameters are set to be 
Nw = 80, Nac = 400 and Vg = 0. 

We first compare the total transmission TLR{e) with 
and without the introduction of the cutoff energy. In 
Fig. [5]^a), we plot TLR{e) as function of the energy e 
with different ac-field strengths Vac for the field fre- 
quency f2 — O.Oltg in the parallel configuration of the 
lead magnetizationsjSS. Here the solid curves represent the 
results without a cutoff, while the dash ones represent the 
results from the calculations by replacing Eq. (jlip with 

g''{er,n,m) = ^ J„-jv(-^)^i-Jv(-^)g5(£r + A^^^) 



N 



X eiE,^t - \er + Nn\) 



(A.l) 



with the cutoff energy E^^t — O.OStg. Our results show 
that after introducing the cutoff energy in such a way, 
the field-free transmission is the same as the one without 
a cutoff when |e| < E^ut, but becomes zero for \£\ > E^^it- 
The difference between the transmissions with and 
without the cutoff energy becomes more pronounced in 
the presence of the ac field. In particular, some pro- 
nounced peaks appear in the transmissions with the cut- 
off energy as indicated by the green and blue arrows, 
which should be absent if calculated properly. This 
phenomenon can be understood via the Tien-Gordon 
theoryfiSi^i^iS which is valid for a weak energy depen- 
dence of the self-energy from the leads. This theory gives 
the formula 



Toie + mQ), 



(A.2) 



where To(e) Tr |f i(£)G5(£)f flXe)G'g(e)| is the field- 
free transmission. The results from this formula are ver- 
ified to coincide with those from our exact calculations. 
This theory shows that the field-applied transmission is 



8 



'cut ^cut 




-0.1 -0.05 0.05 0.1 

E (t,) 



6 



5 



3 



(b) 































0.01 0.02 0.03 0.04 0.05 



FIG. 5: (Color online) Armchair GNR connected with ferro- 
magnetic leads, (a) Total transmission Tlr{s) as function of 
the energy e with different ac-field strength Vac in the parallel 
configuration of the lead magnetizations, (b) Time-averaged 
currents in the parallel and antiparallel configurations as func- 
tion of the ac-field strength Vac- ^ = O.Oltg in these figures. 
The dashed curves represent the results with the cutoff en- 
ergy _Bcut = O.OStg. The blue and green arrows in (a) indicate 
the positions of the pronounced peaks in the blue and green 
dashed curves, respectively. 



just the weighted average of the field-free transmissions 



corresponding to various sidebands. After the cutoff en- 
ergy is introduced, To(£ -I- mfi) in Eq. (|A.2p becomes 
zero for |e -I- milj > -Ecut- Thus, the corresponding side- 
bands cannot contribute to the field-applied transmis- 
sion. When |e| is small, the influence of the cutoff energy 
is still weak, hence the transmission tends to increase 
with increasing |e| just as the one without a cutoff. Nev- 
ertheless, at large |e|, the influence of the cutoff energy 
becomes very important. In this case, the transmission 
tends to decrease with |e| due to the absence of the con- 
tribution of the sidebands satisfying \e + mQ,\ > Ecut- 
As a result, the pronounced peaks in the transmission 
are formed as indicated by the blue arrows. Moreover, 
when the ac field is strong enough, the decreasing trend 
from the cutoff energy is always dominant. Thus, the 
pronounced peak moves to the Dirac point as indicated 
by the green arrow. Based on the above discussions, one 
can conclude that this kind of peaks, which also appear 
in Fig. 2(a) in Ref . [ssl, are just the artificial results caused 
by the introduction of the cutoff energy. 

Similar problems can be found in the field-strength de- 
pendence of the time-averaged current. We plot the time- 
averaged currents I{9) in the parallel {9 = 0) and an- 
tiparallel {0 = tt) configurations against the field strength 
Vac for the Fermi energies Ep = O.OStg and Ep = in 
Fig. IHb) . It is shown that without a cutoff energy, the 
currents in both configurations increase monotonically 
with increasing Vac- This feature is from the increas- 
ing weight of the sidebands outside the energy regime 
between the Fermi energies of two leads, where the field- 
free transmissions are larger than those inside this energy 
regime as shown in Fig. [51[a). However, after the intro- 
duction of the cutoff energy, the behaviour becomes very 
different: the currents first increase and then decrease 
with increasing Vac- The decrease at high Vac is from 
the increasing weight of the sideband beyond the cutoff 
energy, where field-free transmission is zero. This means 
that this decrease, which also appears in Fig. 1(b) in 
Ref. [sl, is an artificial result as well. 
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